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^ , Abstract 

G ■ 

There exist several proposals to use Weizsacker- Williams photons gen- 
erated by ultrarelativistic heavy ions to produce exotic particles in 77 
^ ■ fusion reactions. To estimate the background conditions for such reactions 

we analyze various mechanisms of meson production in very peripheral col- 
lisions of ultrarelativistic heavy ions at RHIC and LHC energies. Besides 
77 fusion they include also electromagnetic 7A interactions and strong 
nucleon-nucleon interactions in grazing AA collisions. All these processes 
are characterised by low multiplicities of produced particles. ^A and AA 
events are simulated by corresponding Monte Carlo codes, RELDIS and 
FRITIOF. In each of these processes a certain fraction of pions is pro- 
duced close to the mid-rapidity region that gives a background for the 77 
events. The possibility of selecting mesons produced in 77 fusion events 
via different pt cut procedures is demonstrated. 
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1 Introduction 



According to the impact parameter, b, different phenomena take place in collisions 
of ultrarelativistic heavy ions. They can be divided into the following three 
categories. 

Central collisions (b m 0), i.e. the collisions with nearly full nuclear overlap 
fall into the first category. Such collisions provide conditions for the creation of 
very hot and dense nuclear matter. It is the aim of future experimental programs 
at Large Hadron Collider (LHC) at CERN |IJ and Relativistic Heavy Ion Collider 
(RHIC) at Brookhaven National Laboratory to study a possible phase tran- 
sition of nuclear and hadronic matter into the so-called quark-gluon plasma at 
high energy densities. As it is expected, the quark-gluon plasma can be created 
under such extreme conditions which are similar to those existed in the Early 
Universe soon after the Big Bang. 

The second category contains collisions with partial overlap of nuclei (b < 
Ri + R2, Ri and R2 are the nuclear radii). In such collisions the residual spec- 
tators remain relatively cold. Short-range interaction via strong nuclear forces 
is restricted mainly to the participant zone. In the whole set of minimum-bias 
events the number of peripheral nuclear collisions is significant due to the geo- 
metrical factor 2nb. The general picture of ultrarelativistic heavy-ion collisions 
at LHC and RHIC will be incomplete without considering such peripheral col- 
lisions. Non-central heavy-ion collisions are considered as a place to look for 
the disoriented chiral condensates || and elliptic flow Combining Hanbury- 
Brown-Twiss method with a determination of the reaction plane || one is able to 
study the size, deformation and opacities of the particle emission sources. With 
the aim to study the central collisions, one has to know the exact properties of 
complementary peripheral collisions. A proper rejection of these "background" 
collision events will be crucial in the extraction of events in which the quark-gluon 
plasma may be created. 

In the collisions of the third category the impact parameter exceeds the sum 
of nuclear radii (b > Ri + R2) . Therefore, there is no overlap of nuclear densi- 
ties, but nevertheless, one or both nuclei may be disintegrated by the long-range 
electromagnetic forces. This process of Electromagnetic Dissociation (ED) is 
a well-known phenomenon |6|, [F], |8|]. The interaction can be treated in terms of 
equivalent photons representing the Lorentz-boosted Coulomb field of heavy ions. 
At ultrarelativistic energies the ED cross section exceeds considerably the pure 
nuclear cross section for heavy colliding nuclei. This fact was confirmed recently 
in experiments ||. The electromagnetic collision events are less violent then the 
collisions with nuclear interactions. Thus the average particle multiplicities are 
essentially lower |Tt], |TJ] and the main part of nucleons and mesons is produced in 
the regions of projectile and target fragmentation, very far from the mid-rapidity 
region. 

Besides the action of virtual (equivalent) photons on colliding nuclei, the pho- 
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tons from two nuclei can collide (fuse) and produce various secondary particles. 
Photon-photon (77) physics may be investigated in such collisions 0, || [12|]. The 
idea to produce exotic particles via 77 fusion in heavy-ion colliders has been put 
forward more then 10 years ago in Ref. ||13|| . Many authors have further devel- 
oped this idea since that time. The production of different particles from /r* 1 , r ± 
leptons to Higgs bosons and supersymmetric particles has been considered. The 
full list of references can be found in the recent reviews [|7|. || . 

As one can see, there exist several mechanisms of particle production in pe- 
ripheral collisions of ultrarelativistic heavy ions: 77 fusion, 7 A interaction and 
grazing nuclear collisions (see Figs. In the present paper we study the 

production of and ir° mesons via these mechanisms. We discuss the general 
features of pion production and calculate the contributions of different mecha- 
nisms. Prior to investigating exotic particle production in 77 collisions a certain 
"calibration" is necessary for the theoretical methods and experimental tech- 
niques. 

Similar investigations have been performed specially for STAR detector at 
RHIC [|T4| , p~5|1 and corresponding acceptance cuts were applied from the begin- 
ning. The CMS detector which will be installed at LHC [16| can also be used for 
studying the meson production in peripheral collisions. Even having in mind a 
plan to use another experimental set-up for such studies one can estimate first 
77 signal to background ratio without any acceptance cuts. Such comparison of 
different mechanisms provides a guide-line which is free of limitations and restric- 
tions of existing experimental facilities. This is useful for possible extensions and 
updates of existing detectors. 

In the present paper we pay the main attention to the inclusive cross sections 
of meson production by different mechanisms. This study is complementary to a 
very recent one |T7j, which deals with the exclusive meson production channels. 
As we expect, in order to select the rare events with the single produced vector 
meson |L7j], one has to reject the background due to meson production (tt + , n~ , 
7T°) via jA process. This background may be estimated in the framework of the 
models used in the present paper. 

A meson production event may be followed by a partial disintegration of a 
nucleus, i.e. a branchy process may take place. We study such reactions in 
addition to the coherent meson production process considered in Ref. |17| when 
the colliding nuclei remain intact. Different characteristics of pion production 
by the equivalent photons are calculated by the RELDIS code |T]J based on 
the extended model of photonuclear reactions |18| . The FRITIOF Monte Carlo 



event generator, version 7.1 |19| is used to study the properties of grazing nuclear 
collisions. 
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2 Distant electromagnetic collisions 



2.1 Equivalent photon spectra 

Let us consider the point-like charge Ze moving with velocity v at the impact 
parameter b. In the Weizsacker- Williams approximation || [7], || the spectrum of 
equivalent photons is given by: 

% 6 ) = ^(* + ^«) W 

where x = cub/'jv, (3 = v/c and 7 = (1 — /3 2 ) -1 / 2 is the Lorentz factor of the 
moving charge. Kq and K\ are the modified Bessel functions of zero and first 
order. 

The total number of photons with energies u)\ and uo 2 colliding at the point P 
is obtained by the integration of the equivalent photon spectra over the distances 
61 and b 2 @ (see Fig. §: 

poo rTm 
IR 2 



roc roo r27V 

F(lu 1 ,lu 2 ) = 2tt bxdbt b 2 db 2 rf0Ar(cu 1 ,6 1 )AT(cu 2 ,6 2 )e(5 2 ), (2) 
JRi Jr 2 jo 



where R\ and R 2 are the nuclear radii, is the step function and B 2 = b\ + b\ — 
2b 1 b 2 cos<p-(R 1 + R 2 ) 2 . 

The 77 luminosity calculated in the double equivalent photon approxima- 
tion |gg is: 

= —F(—e y ™e-y) 

where the energy squared W 2 = Auj\uj 2 and the rapidity y = 1/2 ln(cji/u;2) of the 
77 system were introduced. 

The transverse momentum of produced single meson turns out to be small pt < 
1/-Ri,2 H 0- Different nuclear charge formfactors may be used in calculations || 



|T2| , but in any approximation the p t values for a produced meson turns out to 
be less then 30 MeV/c in PbPb collisions. In the following we shell see how this 
feature may be used to disentangle 77 events from other processes. 

2.2 Production of a single meson in 77 fusion 

The production of a single meson in the 77 fusion is the simplest process in our 
consideration. The calculation technique within the framework of the Weizsacker- 
Williams formalism for such a process is well-known [§, ^] . 

The cross section to produce a meson with mass Mr is given by: 

0"= /— / — F(u 1 ,u 2 )a iy -,M R , (3) 

J Lxl\ J UJ 2 



4 



where the number of colliding photons, F(uji,uj2), is taken from Eq. (0). Taking 
= 1 one can simplify the calculation, since in this case the integral of Eq. (Q) 
may be reduced to the product of the Weizsacker- Williams spectra, N(uj), inte- 
grated over the impact parameter || [7|, ||. For light mesons, Mr <C 7/(i?i + R2), 
the resulting integral is changed by several percent only. This simplification may 
be not appropriate for heavy mesons, since they can only be produced by a pair 
of photons from the high-energy tail of the equivalent photon spectrum. Such 
mesons are produced mainly in close collisions, where b ~ (Ri + -R2) and the 
region of nuclear overlap cannot be neglected. 

The cross section of elementary process 77 — > Mr may be calculated |7|, §, ^0 

as: 

a^M R = 8n 2 (2J R + 1)T Mr ^S(W 2 - M 2 R )/M R , (4) 

where Jr, Mr and 1^/^77 are the spin, mass and two-photon decay width of 
the meson R, W is the cm. energy squared of the colliding photons. 

The corresponding differential cross section for producing a meson with mass 
Mr is given by 

^ = 8^(2^ + 1)F M _ 77 F(^V, fe-)/Ml (5) 

The exclusive cross sections for producing a singe meson were calculated by 
many authors. In Ref. J7| the equivalent photon cross section has been derived 
directly from the first QED principles. Three types of nuclear formfactors were 
used in the calculations corresponding to a homogeneously charged sphere, 
a Gaussian-shaped and a point-like charge distributions. In this respect the ap- 
proach of Ref. |7|] is different from the more phenomenological method of Ref. PU| 
which we basically follow. Moreover, even following the authors of Ref. |20[ one 
can use different values for the nuclear radii, R\^- The real nuclear charge distri- 
bution with a diffuse boundary should be approximated by a distribution with a 
sharp boundary that leads to some uncertainties. A straightforward comparison 
of the results of different authors [|7], ||, ^] seems to be difficult due to different 
values of Mr, r 77 _ > jvf H and R\^ used in the different papers. 

With the aim to understand the sensitivity of numerical results to the choice 
of parameters we repeated the calculation of the cross sections for single meson 
production in ultrarelativistic heavy-ion collisions, see Tabs. |I|-|3|. We used Mr 
and rM H ^77 values according to the corresponding papers 0, ||, pOf| . All the cross 
sections were calculated for the nuclear radii R\^ = 1.2/U. 

As one can see from Tabs. the different approaches give quite similar 
results for the lightest mesons 7r°, 77. On the contrary, as it was expected, the 
difference is noticeable for heavy mesons like rjb, r] c . Nevertheless this situation 
seems to be acceptable due to the following reasons. First, in the present paper 
we will consider mainly the production of the lightest it mesons. Second, be- 
cause of uncertainties in Particle Properties Data on Mr and Tm^^ for heavy 
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mesons [[35[][] it is difficult to predict with confidence the exclusive cross section 
values for such mesons. Therefore, the uncertainties of the method itself become 
less important for heavy mesons. 



2.3 Inclusive cross section of rr production in the DRP 
model 

Pions can be produced not only in the direct process 77 — > n°. The 77 fusion 
produces several unstable heavy mesons such as r],f ,... and decay products 
of these mesons may contain 7r°'s. As one can see in the following, such two- 
step processes, 77 — > R — > ir°X, with the intermediate meson R give a sizable 
contribution to the inclusive tt° production cross section. In addition to the 
direct process a reasonable estimation for the inclusive cross section 77 — ► 7r°X 
should contain the sum of dominant contributions from the intermediate meson 
resonances which contain at least one 7r° meson in the final state. We call this 
approach a Dominant Resonance Production (DRP) model. 

The exclusive differential cross section da^jdy to produce a single meson R 
with mass Mr is given by Eq. ([5]). In the DRP approximation the inclusive cross 
section of the 7r° production through the decay of intermediate resonances is 

^V) = E ^ b r( r - -°) »2V). (6) 

where B^(R — > 7r°) is the branching ratio for decay of the resonance R to 
the channel k which contains at least one 7r°. n^(7r°) is the number of 7r°'s in 
the corresponding channel k. If the first step decay products contain another 
resonance, a similar expression may be written in turn for the second step decays 
producing 7r°. Introducing the value 

B° R ut = Y, 4 fc) (*-TV2V) (7) 
k 

one can rewrite Eq. @ in the following way: 

d<Ti nc l / n\ R uout /n\ 

■*- ( ' ) = §4T B «- (8) 

Particle Data information relevant to the calculation of Oi nc i is given in Tab. [|. 
We selected the resonances with relatively large widths Yu R -yrt- ^ e vames °f ^tot 
and = YM R -*rr/Ytat from review |^5| which are necessary to calculate rjvf H -> 77 
are given for completeness. On the contrary to Tabs. the decay probabilities 
and meson widths in Tab. ^ were taken according to the most recent version of 

1 This is particularly true for scalar mesons /o(975) and /o(1275) with large decay widths 
causing a strong overlap of individual resonances, see details in Ref. [p5| . 
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Review of Particle Physics [p^] . We assumed that the decays /o(980) — > tttt and 
/o(1370) — > mi take place with 100% probability since there is no quantitative 
information in Ref. []35] on other channels. Isospin conservation relations provide 
the probability of 1/3 for charge state in such tttt channels. 

The integrated inclusive cross section of 7r° production in PbPb collisions 
(Ri = i?2 = 7.75 fm) at LHC energies is found to be cxj nc ;(7r ) = 106 mb in DRP 
model, while the exclusive one is only 36.25 mb. As one can see, the contribution 
to the 7T° production via intermediate resonances turns out to be essential. The 
corresponding rapidity distributions, da/dy, will be discussed in Sec. §. 

2.4 Pion production in collisions 

The meson production induced by equivalent photons in electromagnetic colli- 
sions of ultrarelativistic heavy ions (Fig. is a poorly explored phenomenon. To 
the best of our knowledge, the first calculations of the total rate of pion produc- 
tion in electromagnetic collisions were made in Ref. [22|. The first experimental 



evidence of electromagnetic dissociation accompanied by the pion production was 
found in Ref. [23[] for 200 AGeV 16 ions. Due to a small number of pion pro- 



duction events detected in nuclear emulsion, the absolute cross section for such 
dissociation channel was not determined. To date only the existence of pion 
production in 'jA collisions is demonstrated by the experiment |2"3]| . 

Let us consider the absorption by a nucleus of an equivalent photon leading to 
the pion production (Fig. ^j. The nucleus may absorb one or more virtual photons 



during a collision. We follow Llope and Braun-Munzinger |2^] in description of 
such multiple absorption processes. The double differential cross section of the 
pion production via the single photon absorption may be written first in the 
projectile rest frame. In such a frame the nucleus is at rest prior to absorption: 

1 1 -N^(u 1 )a A2 (u 1 )—-—(u 1 ), (9) 



Ptdp t dy J uji Ptdptdy 

where oa 2 is the total photoabsorption cross section for the nucleus A 2 and 
d 2 W / ptdptdy is the double differential distribution of pions produced by a photon 
with energy oj\. Spectral function is given by the following expression: 

N {1 \u 1 ) = 27r J bdbe'^Niuub), (10) 

where N(ui,b) is defined by Eq. (JJ) and m(b) is the mean number of photons 
absorbed by the nucleus A 2 in a collision at impact parameter b: 

m(b)= [ N(co,b)a A2 (u)—. (11) 
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Analogously, for the second order process with a pair of photons absorbed by the 
nucleus A 2 the double differential cross section is: 



(12) 



Ptdptdy J J loi lu 2 Ptdptdy p t dp t dy 

with the corresponding double photon spectral function: 

00 

N^{uj x ,uj 2 ) =tt f bdbe- m(b) N{uj 1 ,b)N{uj 2 ,b). (13) 



In the above expressions b m i n is the minimal value of the impact parameter 
which corresponds to the onset of nuclear interaction. In pion production cal- 
culations the integration over the photon energy starts from LO m i n Pd 140 MeV 
which corresponds to the pion emission threshold. Finally, with the correspond- 
ing Lorentz boost from the nucleus rest frame to the laboratory system one can 
obtain the double differential distribution d 2 a/p t dp t dy for produced pions suit- 
able for measurements in experiments. Further details of our approach may be 



found in Refs. llOL 11 



As it was shown in Ref. [OJ|, the contribution to pion production from the 
double photon absorption is less then 10% for heavy ions with 7 ^> 100. Since the 
contribution from the third and fourth order processes are expected to be even 
lower, these processes can be safely disregarded for the ultrarelativistic energies. 

Because of the coherent action of all the charges in the nucleus, the virtu- 
ally of the emitted photon, Q 2 = —q 2 , is restricted. Such photons are almost 
real, Q 2 < 1/R 2 , where R is the nuclear radius. Therefore, photonuclear data 
obtained in experiments with monoenergetic photons may be used, in principle, 
as an input for Weizsacker- Williams calculations of pion production in jA colli- 
sions. Since the spectrum of equivalent photons covers a very wide range of the 
photon energies, one needs the double differential distribution d 2 W (a>) / p t dp t dy 
for photon energies uj starting from the pion emission threshold and up to several 
tens or even hundreds of GeV. 

In the region of interest, i.e. at u > 140 MeV, the photon de Broglie wave- 
length is comparable or even smaller than the nucleon radius. A photon interacts, 
mainly, with individual intranuclear nucleons. Experimental data on single pion 
photoproduction on the nucleon are accumulated in compilations of Refs. [^5|, . 
The latter compilation contains also the experimental data on photoproduction 
of baryon, B*, and meson M* resonances: 'yN — > nB* and 7^ — > tcM* as well 
as on some channels of multiple pion production: jN — ► iirN, 2 < i < 8. Due 
to a long mean free path an equivalent photon may be absorbed deeply in the 
nuclear interior. The mesons produced in a photonucleon reaction interact with 
other intranuclear nucleons inducing different reactions in the nuclear medium. 
In other words, the process shown in Fig. |2| is followed by the final state interac- 
tion of produced hadrons with the residual nucleus A 2 . Beside the pions several 
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nucleons may be emitted and the residual nucleus may receive a high excitation 
energy. 

Since the data on pion photoproduction on nuclei exist only for limited do- 
mains of to and pion kinematical variables, a theoretical model should be used 
in the cases when the data are not available. A suitable tool for describing such 
multi-step photonuclear reactions is the Intranuclear Cascade Model, which is 
well-known for many years ||27|1 . 



As shown in Ref . |T8[ , the extended Intranuclear Cascade Model of photonu- 
clear reactions describes reasonably well available data on meson production and 
nucleon emission obtained in last two decades with intermediate energy quasi- 
monochromatic photons. With this in mind one can use the INC model for the 
Weizsacker- Williams calculations, as it was demonstrated in Refs. [1^, DJ. For- 



mally we use 1 TeV as the upper limit for the equivalent photon energy in the 
7 A process. Our model was not initially designed for such high energies when 
events with a very high hadron multiplicity become possible and many other 
channels, like baryon-antibaryon photoproduction, may be open. Nonetheless 
one can safely use the model to investigate the jA processes with low multi- 
plicity which mimic processes from 77 fusion. The former processes take place 
mainly at u < 10 GeV, where our model has been verified in detail. 

2.4.1 Simulation of 'jN interaction 

The channels of the hadron photoproduction on a nucleon which are taken into 
account in the model are listed in Tab. |5|. To describe the two-body photo- 
production channels we have basically used the Monte Carlo event generator of 
Corvisiero et al. | 2q| . The contribution from the two-body channel 7iV — > ttN 



dominates up to uj ~ 0.5 GeV. Approximations of the total and differential cross 
sections for such channels based on the model of Walker and Metcalf |29| were 



used up to 2 GeV. The excitation of six baryon resonances was taken into account. 
The contributions from A(1232), iV*(1520) and iV*(1680) are the most impor- 
tant. The presence of these resonances explains a well-known resonant structure 
in the total 7^ cross section at to < 1.2 GeV. 

The channels 7iV — > 2tcN and 7iV — > 3ttN play a major role at 0.5 < uj < 
2 GeV. These channels include also the resonant contributions from tcA, r)N, 
pN and ujN channels (Tab. ||). Although the presence of these contributions 
in the total 'yN cross sections is difficult to trace, the angular distributions of 
these channels have some specific features. Several examples may be given. The 
7P — * t]p channel has three important contributions from 5 , n(1535), 5*11(1700) 
and Pn(1750) states. Due to the dominance of 511(1535) the angular distribution 
of the process is not far from isotropy. On the contrary, 7p — > p°p process 
has a prominent forward peak, since the non-resonant diffractive contribution 
dominates. 

Many channels in the multi-pion reactions, 7iV — > iirN, (2 < i < 8), were not 
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suitable for measurements in spark, bubble or streamer chamber experiments due 
to the presence of several neutral particles in the final state. However, one can 
reconstruct the integral cross sections of undetectable channels by applying iso- 
topic relations to the measured cross sections of channels with charged particles 
which can be found in compilation of Ref. [2B|. A phenomenological statisti- 
cal approach for the exclusive description of the elementary 7 A" interaction was 
applied to multiple pion production channels in Ref. ||18|| . There an isospin sta- 
tistical model was used to connect unknown cross sections with measured ones. 
Other details of the method used for simulating the 7 A" interaction may be found 



elsewhere, see Ref. 18 



Since a huge number of multiple pion production channels is open at to > 2 
GeV, the statistical description may be the only way to estimate the cross sections 
of such channels. Recently another kind of statistical assumptions was used in 
Monte Carlo event generator Sophia for simulating photohadronic processes in 
astrophysics [KO . 



2.4.2 Secondary processes 

According to our model the fast hadrons produced in a primary 7 A" interaction 
initiate a cascade of successive quasi-free hadron-nucleon collisions inside the 
nucleus. The following elementary processes were taken into account: 

vrN -> vrN , n(NN) NN , vrN -> vrvrA^ ; vrN -> (i + l)nN , 

NN -»• NN , NN -> irNN , NN -> inNN , (i > 2) ; 

7]N -> 7]N , 7]N # nN , rjN -> rnrN , rj(NN) -> NN , rf(NN) -> ttA^AT , 

cuA^ -> ujN , wiV # vrA^ , uoN -> ttttA^ , w(iViV) -> iViV , w(JVJV) -> TriViV . 

The empirical approximations for the measured integral and differential cross 
sections of A^V and nN interactions as well as the phenomenological estimations 
for the total and partial cross sections of r)N and ujN interactions were used in 
the calculation llSfl. 



The described Monte Carlo model is implemented into the RELDIS code [11 



which is especially designed for calculating the electromagnetic dissociation pro- 
cesses in ultrarelativistic heavy ion collisions. 



2.4.3 Characteristics of produced pions 

The average number of pions, (n^), produced in a jA process is quite small. 
The values of (nN for each pion charge are given in Tab. ^| for AuAu and PbPb 
collisions at RHIC and LHC energies, respectively. Even with the inclusion of 
multiple pion production channels (7 AT — ► iirN, 2 < i < 8), the average numbers 
of pions of each charge are quite low (w 1). This can be explained by two 
reasons. First, in a great part of the electromagnetic dissociation events no pions 
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are produced at all, due to the dominance of soft photons (u < 140 MeV) in the 
equivalent photon spectrum, Eqs. (|T0D , ([13|) . Second, the most probable processes 
of pion production are 7iV — > 7riV and 'jN — > 71 A, which produce one or two 
pions only. 

As one can see from Tab. |^, in the 7A process the neutral pions are produced 
more abundantly as compared with tt + and tt~ mesons. A difference is noticeable 
also in the differential distributions, da/dy, shown in Fig. ||] for each pion charge. 
The main deviations in the pion yields appear at rapidities close to the beam 
rapidity. This fact has a simple explanation. As already mentioned above, the 
process 'jN — > ttN dominates in the pion production by equivalent photons and 
corresponding pions are close to the beam rapidity. This reaction can proceed 
on a proton: 7£> — > ir + n, 7p — > 7r°p and a neutron: — > vr~p, — > 7T°n. 
Thus, single neutral pions can be produced both on the proton and the neutron, 
while 7r + on the neutron and 7T~ on the proton only. The total cross sections of 
these four channels are close to each other. This gives the rate of 7r° production 
approximately twice as large as the rate of 7r + or 7r~ for a light nucleus with 
the equal numbers of protons and neutrons. This feature of the single pion 
photoproduction is confirmed by the measurements of inclusive cross sections 



and the calculations made for carbon nucleus (see Fig. 15 of Ref. [|T| 



Since a neutron excess exists in heavy nuclei like Pb or Au, in this case the 
7r~ production will be enhanced compared to tt + , as it is shown in Fig. [5[ The 
final state interaction affects the absolute yields of pions of different charges, but 
does not essentially change the ratio between them. 

The detection of pions produced close to the beam rapidity is a complicated 
experimental task. Special zero-degree detectors with a proper identification of 
particle mass and charge are necessary for this purpose. However if a forward 
detector located after a steering magnet is suitable for determination of charges 
of nuclear fragments, one may exploit it to detect the pion emission process in- 
directly. If 7r~ is produced in the reaction jn — > n~p, it may be emitted while 
the recoil proton may be captured by the residual nucleus. The deexcitation of 
such nucleus may take place mainly via neutron emission. It means that the 
initial charge of the nucleus will be increased by one unit. Charge-to- mass ratio 
will be changed for such ion and it will be separated from the beam. Depend- 
ing on the value of heavy-ion energy RELDIS code predicts the cross section of 
such electromagnetic "charge pick-up" channels at the level of 10-100 mb. This 
process competes with the nucleon pick-up process via strong interaction. Nev- 
ertheless as it is known from Ref. |31|, the cross section of the latter process 
drops noticeably with increasing beam energy. On the contrary, the cross section 
of electromagnetic charge pick-up increases gradually with the beam energy and 
may essentially exceed the value of cross section for the nucleon pick-up due to 
the strong interaction. 

The photonuclear reactions (7, 7r~xn), x = — 9 induced by bremsstrahlung 
photons were discovered many years ago. We refer the reader to a recent paper 
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where such reactions were studied by radiochemical methods |32|. It would be 
interesting to study the same type of reactions induced by equivalent photons in 
ultrarelativistic heavy ion collisions. 

As mentioned above, the main part of pions is produced close to the beam 
rapidity A small fraction of pions is produced by a very-high energy photons 
with uj ^> 10 GeV. Some of the pions produced in such interactions may receive 
the momentum high enough to populate the central rapidity region. As shown 
in Fig. |5], there are pions with \y\ < 5 which may be confused with those from 
77 process. Therefore some selection criteria are necessary for 77 processes (see 
Sec. W). 



3 Very peripheral nuclear collisions 

Very peripheral (grazing) nuclear collisions with the participation of strong nu- 
clear forces can be misidentified as the electromagnetic interaction events. Both 
of these types of interactions contribute to the events with low multiplicity of 
particles. To study the properties of strong interaction events a Monte Carlo 



event generator FRITIOF, version 7.1 [19] is used 



3.1 Basic processes considered by FRITIOF model 

Let us recall the main statements of the extended FRITIOF model |T^| which 
is valid up to the TeV region. FRITIOF is a Monte Carlo model for hadron- 
hadron, hadron-nucleus and nucleus-nucleus collisions. The basic idea of the 
model is in a simple picture that a hadron behaves like a relativistic string with a 
confined color field. This field is equivalent to that of a chain of dipoles lined up 
along the axis line. The dipole links act as partons. During the soft interaction 
many small transverse momenta are exchanged between the dipole links and two 
longitudinally excited string states result from the collisions. A disturbance of 
the color field will in general initialise the gluonic radiation according to QCD. 
The final state particles are obtained by fragmentating the string states like in 
the e + e~ annihilation. 

The large p t process can be treated by using QCD directly. The hard inter- 
action effects, which are considered as the Rutherford parton scattering, become 
important in the TeV range of cm. energies and at large p t > 1 GeV/c. The re- 
sults of the model are in good agreement with experimental hadron-hadron data 
up to the highest energies currently available. 

Nucleus-nucleus collisions are regarded in the FRITIOF model as incoherent 
collisions between nucleons of colliding nuclei. FRITIOF does not take into ac- 
count collective (coherence) effects when two nuclei interact as a whole. Thus a 
nucleon from the projectile interacts independently with the encountered target 
nucleons as it passes through the nucleus. Each of these subcollisions can be 
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treated as a usual hadron-hadron collision. On the time scale of the collision 
process, the exited nucleon does not fragment inside the nucleus, so there are no 
intranuclear cascades. This assumption is reasonable at high energy since the 
time scale associated with fragmentation is much longer than the flight time of 
excited nucleons through the nucleus. 

3.2 Nuclear density distributions and total cross section 
in PbPb collisions 

One of important aspects of the collision is the nuclear geometry. It is assumed 
that the projectile passes through the target nucleus on a straight line trajectory. 
The nucleons are distributed inside the nucleus according to the nuclear density 
distribution p(r). Wood-Saxon density function for heavy nuclei was used in our 
calculations: 

p(r) = ^ (14) 

where r Q = 1.16(1 — 1.16A -2 / 3 ) is the radial scale parameter, C is the diffuseness 
parameter, it is slightly A-dependent and its values used in the program ranging 
from 0.47 to 0.55 fm, and p Q is a normalisation constant. Then for lead nucleus 
the radius defined at the half of the normal nuclear density is equal to R = 
ToA 1 / 3 = 6.63 fm, while the diffuseness parameter used by FRITIOF is equal to 
C = 0.545 fm. 

The total cross section of strong lead-lead interaction at LHC collider was 
estimated to be 7164 mb, according to a simple geometrical formula used in 



Ref. [Q. This estimation does not take into account the energy dependence of 
the total nucleus-nucleus cross section due to the increase of the total nucleon- 
nucleon cross section. Basing on extrapolations of the measured pp cross sections 



to LHC energies one can obtain a higher value of about 10 barn, see Ref. [p4 
For this value of the total nucleus-nucleus cross section one has to rescale the 
plots given in the present paper by 30% for grazing nuclear collisions. 



3.3 Impact parameter distribution in nuclear collisions 

Special efforts should be undertaken to find proper ranges of the impact pa- 
rameter b which divide nucleus-nucleus collisions into different categories. These 
categories were listed in the introduction, most important for us are very periph- 
eral nuclear collisions and electromagnetic interactions, subsequently subdivided 
into 77 and 7 A events. 

The uncertainties associated with the diffuseness of the nuclear density dis- 
tribution affect grately the resulting ranges. We used the FRITIOF model to 
define these ranges. The impact parameter distribution of events simulated by 
the FRITIOF code is given in Fig. |6|. The plotted quantity dN aa(P) / db / N aa 
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is a relative number of nuclear collisions with given b. As it is shown in the 
figure, the events with a low multiplicity of pions correspond to collisions with 
12 < b < 20 fm. It is the domain of b which defines the grazing nuclear collisions 
with participation of strong nuclear forces. It should be stressed however, that 
this definition is essentially model dependent and numerical results may be dif- 
ferent if one uses another code with different treatment of the diffuse boundary 
of the nucleus and another values of the elementary NN cross sections. 

We see that the value which defines the onset of electromagnetic collision 
for PbPb ions, b m i n = 15.5 fm, used in the paper |10| does not contradict the 
results of the FRITIOF calculations given in Fig. 10. Following our previous 



papers |T0| , [TT[ , for consistency we define the end of nuclear overlap in lead-lead 



collisions at b m i n = 15.5 fm. This choice is also supported by the calculation of 
Ref. [36 where the average number of interacting nucleons, (N in ), was evaluated 
as a function of b for the case of PbPb collisions. As it was found J5BJ, (N in ) « 1 
at b = 15.5 fm. According to the Wood Saxon density function for Pb nucleus it 
corresponds to the overlap of those boundary regions of the colliding nuclei which 
have the nuclear densities below 0.1 of the value at the center of the nucleus. The 
total number of such events in the FRITIOF simulation turns out to be 24%. 

Fig. |7| demonstrates the distribution of pion multiplicities for peripheral PbPb 
collisions (b > 15.5 fm) at LHC energies simulated by the FRITIOF code. The 
most probable is the production of four pions and the distribution is very broad. 
This conclusion is valid for all the charge states of pions. The average pion 
multiplicities are given in Tab. where one can see that the multiplicities of 
pions produced in the above-defined grazing AA collisions are generally greater 
than those for 7 A process. 



4 Comparison of different mechanisms of pion 
production 

Since most of the detectors used thus far in heavy-ion experiments have very 
restricted acceptance on rapidity, it is important to investigate the rapidity de- 
pendence of pion production cross section. Rapidity distributions of neutral pions 
produced in peripheral PbPb collisions at LHC energies are shown in Fig. [8[ 

Let us consider first the distributions without any additional selection criteria. 
The distributions from 77 fusion and grazing nuclear collisions calculated by the 
FRITIOF code with b > 15.5 fm are very similar in shape. The distributions of 
da/dy for single and inclusive ir° production in 77 fusion also have maximuma at 
y = 0. But even in this region it is much lower then the contribution from grazing 
nuclear collisions. On the contrary, the rapidity distribution for 7A production 
has peaks at the beam and target rapidities. Nevertheless its contribution at 
midrapidity is comparable to the one from 77 fusion. 
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This above-described picture is unfavourable for experimental detection of 
pions from 77 collisions. Nonetheless one can improve the detection conditions 
by using appropriate transverse momentum cuts. Indeed the single mesons from 
77 fusion have a very small transverse momenta. Selecting events with low values 
of pt one can reject contributions from other processes. Following this way two 
different p t -cut procedures may be proposed. 

According to the first procedure (b-criterion) one should select the events with 
a small total transverse momentum of the meson system \p t {sum) \ < 75 MeV/c. 
The second selection procedure (c-criterion) imposes more severe restrictions on 
the transverse momenta of the produced pions demanding each of them (including 
7T°, 7r + , and 7r~) to be small, \p t \ < 75 MeV/c, in this case. As one can see, b- 
criterion is equivalent to c-criterion in the case of single pion production. 

Let us investigate now to what extent the rapidity distributions are affected 
by p t cuts. The distributions obtained according to b- and c-criteria are also 
given in Fig. |8|. One can see that the b-criterion reduces the contribution from 
grazing nuclear collisions by three orders of magnitude and give some benefits 
for detection of 77 events. If this reduction is not sufficient one can use c- 
criterion which is more efficient. Further suppression of pions from grazing nuclear 
collisions may be obtained in the region —4 < y < 4. With this suppression even 
exclusive 77 — > 7r° channel may be clearly distinguished. 

The difference in b- and c-criteria applied to 7 A events calculated by RELDIS 
code is less prominent. This is explained by the fact that single pion production 
dominates in 7 A collisions at u > 140 MeV while b- and c-criteria are equiva- 
lent for this case. Both of the procedures may be recommended to reduce the 
contribution from 7 A process as it is shown in Fig. KL 

It should be noticed that b-criterion does not affect the inclusive distribution 
for 77 fusion, while c-criterion suppresses the main part of events of 7r° production 
due to heavy meson decays. This may be useful for extraction of the exclusive 
77 — > 7T° process. 

5 Conclusions 

Only several hadrons are produced on average in very peripheral collisions of 
ultrarelativistic heavy ions. Because of this feature such collisions can be con- 
sidered as non-violent events. An exact determination of the impact parameter 
in a collision event is beyond present experimental techniques. However, as it 
is shown in the present paper, selecting the events with very low multiplicity 
of produced mesons one can approximately identify a domain of large impact 
parameters where 77 fusion, 7 A or grazing nuclear AA collisions takes place. 

Calculations show that each of the mechanisms has specific distribution of 
produced pions on transverse momentum, p t , and rapidity, y. One can use these 
features for adjusting detectors to different regions of p t and y with possibility 
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to disentangle pions produced by different mechanisms. One can enhance the 
signal to background ratio for 77 fusion by selecting events with low p t by means 
of two procedures considered in the paper. Our preliminary results confirm the 
importance of p t cut proposed in Ref. fL2| . This procedure rejects the background 



from 7 A and grazing AA collisions by several orders of magnitude. Moreover 
using p t cuts one can avoid a pessimistic conclusion made in Ref. |T7| that the 



77 fusion is indistinguishable from other process since the rapidity distributions 
for different processes are very similar to each other. 
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Table 1: Exclusive cross sections of single meson production in 77 fusion process 
in ultrarelativistic UU and PbPb collisions at colliders. The values were calcu- 
lated for comparison with Ref. [pO|l . The masses and two-photon decay widths 



for pseudoscalar mesons were used according to Ref. 
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Table 2: Exclusive cross sections of single meson production in 77 fusion pro- 
cess in ultrarelativistic PbPb collisions at LHC The values were calculated for 
comparison with Ref. ||. The masses and two-photon decay widths for cc and 
bb mesons were used according to Ref. ||. 
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Table 3: Exclusive cross sections of single meson production in 77 fusion process 
in ultrarelativistic AuAu and PbPb collisions at RHIC and LHC The values were 
calculated for comparison with Ref . . The masses and two-photon decay widths 
for different scalar and pseudoscalar mesons were used according to Ref. 0. 
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Table 4: Resonances, their decay modes and branching ratios (given in brackets) 
which were taken into account in inclusive cross section calculation in DRP model. 
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Table 5: Channels of elementary yN interaction taken into account in the INC 
model. 
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Table 6: Average numbers of tt + ,tt and ir° mesons produced in yA and grazing 
AA collisions 
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Figure 1: Meson production in 77 fusion. 




Figure 2: Meson production in jA collision. 
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Figure 3: Meson production in grazing AA collision. 




Figure 4: Effective 77 luminosity for heavy-ion collisions. Beam direction is 
perpendicular to the picture plane, bi and 6 2 are the distances from the nuclear 
centers to the photon interaction point P. 
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Figure 5: Rapidity and transverse momentum distributions of pions produced in 
'-/A process in PbPb collisions at LHC. The results of RELDIS code are given 
by solid, dashed and dotted histograms for 7T + , tt~ and tt° mesons, respectively. 
Beam rapidities are shown by arrows. 
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Figure 6: Impact parameter distributions for nuclear PbPb collisions at LHC 
simulated by FRITIOF code. The solid histogram corresponds to all events. 
The doted and dashed histograms present the impact parameter distributions for 
events with N n o = 1 and N n o = 2, respectively. 
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Figure 7: Pion multiplicity distributions for grazing nuclear PbPb collisions at 
LHC simulated by FRITIOF code. The events with impact parameter b > 15.5 fm 
were selected. 
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Figure 8: Rapidity distribution of 7r° produced in peripheral PbPb collisions 
at LHC energies. The thick solid and dashed line histograms labeled "1" are 
the inclusive (DRP model) and exclusive cross sections for 7r° production in 77 
fusion. The histograms "2a" , "2b" and "2c" give the results of RELDIS code for 
7 A process. The FRITIOF results for peripheral PbPb collisions with b > 15.5 fm 
are given by the points "3a" , "3b" and "3c" . The distributions with label "a" were 
obtained without p t cuts. The label "b" corresponds to the selection of events 
with the total transverse momentum of meson system |p^ smn )| < 75 MeV/c. The 
label "c" corresponds to the selection of events with \p t \ < 75 MeV/c for each of 
the pions including n°, 7r + , and n~ . 
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